Measurements of the Decay K L — > e 
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The E799-II (KTeV) experiment at Fermilab has collected 83262 Kl -> e+e~7(7) events above 
a background of 79 events. We measure a decay width, normalized to the Kl 
(ty" — > 77, 7T° — > 77, tt% — > e + e~7(7)) decay width, of Y{Kl — > e + e _ 7(7))/r(_K"i, - 
(1.3302 ± 0.0046st a t ± 0.01028^34) x 10 -3 . We also measure parameters of two A"_l7*7 form fac- 
tor models. In the Bergstrom, Masso, and Singer (BMS) parametrization, we find Cok* = 
—0.517 ± 0.030 s tat ± 0.022 sys t. We separately fit for the first parameter of the D'Ambrosio, Isidori 
and Portoles (DIP) model and find a DIP = -1.729 ± 0.043 s t a t ± 0.028 sys t. 



PACS numbers: 13.20.Eb, 13.25.Es, 14.40.Aq 

The rare decay Kl e + e~j offers a direct means 
for studying the dynamics of the Kt,j*j vertex. The 
form factor of this vertex is important for determin- 
ing the long-distance, two photon contribution to the 
Kl — > decay width so that the more interesting 

short-distance contributions can be determined. These 
contributions have been important for extracting a di- 
rect constraint on the real part of the CKM matrix pa- 
rameter Vtd [Jill]- In addition, this short-distance infor- 
mation provides one of the most stringent constraints on 
flavor-changing neutral-current couplings of the Z boson, 
which are beyond the standard model In this letter 
we present a measurement of the KlJ*^/ form factor and 
Kl — * e + e~7 branching ratio using data from the 1997 



run of the E799-II (KTeV) experiment at Fermilab. 

Two parametrizations of the form factor are consid- 
ered. The first model, proposed by Bergstrom, Masso, 
and Singer (BMS), describes the i\L7*7 vertex in terms 
of two types of processes: a Kl transition into a ir°,r),r]' 
state that decays to two photons, and a vector meson 
dominance contribution in which the Kl first decays into 
a K* and a photon followed by a strangeness changing 
vector-vector transition into a p,to,4> state that decays 
into a virtual photon The relative size of these two 
contributions is characterized by a constant parameter, 
Ccyk* ■ The BMS form factor model, given in Eq. [U is a 
function of the Dalitz variable x, which is defined as the 
squared ratio of the e + e~ mass to the Kl mass. 
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The form factor model proposed by DAmbrosio, Isidori, and Portoles (DIP) is a slightly more general model that 
applies to all Kl7*7* vertices 

fDip(xi,x 2 ) = 1 + a DIP ( M 2 + J^2 I "b PDIP , M 2 \ f M 2 \ 
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The variables x\ and xi are the squared ratios of the 
masses of the two virtual photons to the kaon mass. The 
Kl — * e + e~7 decay is only sensitive to aoip since one 
of the photons emerging from the vertex is real. 

The E799 phase of the KTeV experiment focused 800 
GeV protons from the Tevatron at Fermilab onto a BeO 
target. The resulting particles were collimated into two 
parallel neutral beams. The beams passed through a 65 
m vacuum decay region beginning 94 m downstream of 
the target. Immediately following the decay region was 
a charged track spectrometer. It consisted of two up- 
stream drift chambers, a dipole magnet, and two down- 
stream drift chambers. The drift chambers achieved a 
position resolution of about 100 /im, which corresponded 
to a momentum resolution of a(P)/P = 0.38% © 0.016%- 
P(GeV/c). Both the decay pipe and charged spectrome- 
ter were surrounded by lead-scintillator photon veto de- 
tectors. 

A transition radiation detector (TRD) consisted of 
eight planes of polypropylene felt, each followed by a 
multi-wire proportional chamber (MWPC) containing an 
80%-20% mixture of xenon and C0 2 . The TRD provided 
a single track pion rejection of better than 200:1 at 90% 
electron efficiency [5J. 

A Csl electromagnetic calorimeter was just down- 
stream of the TRD. The calorimeter provided an energy 
resolution of a(E)/E = 0.45% © 2%/y/E(GeV). Behind 
the calorimeter was a muon system consisting of planes 
of lead, steel and scintillator. More detailed descriptions 
of the KTeV detector can be found elsewhere @, Q . 

A Kl — > e + e~7 event was observed in the KTeV detec- 
tor as two oppositely charged tracks in the charged spec- 
trometer that each pointed to a cluster in the calorimeter, 
and a third calorimeter cluster corresponding to the pho- 
ton. Each cluster was required to have an energy greater 
than 2.75 GeV, and for electrons the ratio, E/p, of cluster 
energy (measured by the calorimeter) to track momen- 
tum (measured by the spectrometer) was required to lie 
between 0.925 and 1.075 to discriminate against muon 
and pion backgrounds. Several cuts were made to en- 
sure that the decay vertex was inside the fiducial region 
and to remove events near the edge of a detector. A ver- 
tex x 2 cu t required that the decay particles originated 
from a common vertex, and a \ 2 cu t on track match- 
ing at the center of the analysis magnet ensured that the 
tracks were well-reconstructed. The total kaon energy 
was required to be between 40 GeV and 200 GeV, and 
less than 0.15 GeV of energy was allowed in each of the 
photon veto counters. 

Since all of the final state particles in this decay mode 
were reconstructed, it was possible to impose two addi- 
tional kinematic constraints. The first was a cut on the 
square of the component of the reconstructed kaon mo- 
mentum transverse to the original kaon direction (pf). A 
cut on p\ was placed at 500 (MeV/c) 2 to reduce events 
with extra particles not related to the decay and events 




Reconstructed e + e y Mass (GeV/c 2 ) 

FIG. 1: The reconstructed e + e~7 mass distribution before 
and after the TRD cut shows a significant decrease in the 
amount of background while retaining most of the signal. The 
signal region is indictated by the two vertical lines centered 
on the kaon mass. 



with missing decay particles. The other kinematic con- 
straint was the requirement that the reconstructed kaon 
mass lie between 0.475 GeV/c 2 and 0.520 GeV 2 . 

One of the main backgrounds was from Kl — > 7r ± e =F ^ e 
(K e 3) decays with either an accidental or radiated pho- 
ton. In order for these events to mimic a signal, the pion 
had to be misidentified as an electron. This occurred 
when most of its energy was deposited in the calorime- 
ter yielding an E/p value close to one. The K e 3 back- 
ground peaked below the signal in reconstructed e + e~7 
mass, which is consistent with the kinematic limit of 
0.477 GeV/c 2 for an event with a pion misidentified as 
an electron. This background was significantly reduced 
by requiring a TRD pion probability of less than 5% for 
each track. Fig. [TJ shows a plot of the e + e~7 invariant 
mass before and after the TRD cut. 

The other main background was due to Kl — * 77 
events where one of the photons converted into an e + e~ 
pair while exiting the vacuum region. Since these pho- 
ton conversions yielded e + e~ tracks that tended to be 
close together upon reaching the charged spectrometer, 
requiring a separation of 1.5 mm in the x or y view of the 
first drift chamber of the charged spectometer effectively 
removed these events. In a Monte Carlo simulation of 
10 times the expected number of Kl —> 77 conversion 
events, no events survived this cut. 

To measure the branching ratio, the Kl — * e + e~7 
rate was normalized to the rate of Kl — > 7r 7r 7r^, 
(Kl — > 7r°7T 7r with two subsequent ir° — > 77 decays 
and one ir° — > e + e _ 7 decay). To reconstruct this mode, 
two tracks pointing to clusters in the calorimeter and five 
additional photon clusters were required. To determine 
which four photons belonged to the two tt° — > 77 de- 
cays, the decay vertex position was calculated for each 
pair of photons using the separation distance and en- 
ergies of the corresponding clusters, and assuming that 
they came from a ir° . A pairing \ 2 w &s then constructed 
based on the requirement that both 7r°particles had the 
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same vertex position. The photon pairing that yielded 
the smallest % 2 was assumed to be correct. The remain- 
ing photon was assumed to come from the ir° — ► e + e~7 
decay. Photon mispairing occurred a small fraction of 
the time, but this effect was well simulated in the Monte 
Carlo. 

The cuts used to isolate Kl — > ■k®-k°it d events were 
kept as similar as possible to those used for Kl — * e + e~7 
with a few exceptions. The four additional photons in 
the final state significantly improved the total kaon mass 
resolution, so the kaon mass cuts were tightened to 0.485 
GeV/c 2 and 0.510 GeV/c 2 . A loose cut was made on 
the photon pairing \ 2 > an d after matching the remaining 
photon with the e + e~ pair, a cut was made on the e + e~7 
mass to select the region between 0.1275 GeV/c 2 and 
0.1425 GeV/c 2 . 

The only significant background to the decay Kl — > 
7r 7r 7r^, was from Kl 7r 7r 7r° events where one of the 
photons converted to an e + e~ pair. This background 
was removed with the same track separation requirement 
used in the signal mode. 

Monte Carlo simulations were used to determine the 
acceptances of the signal and normalization decay modes. 
Both modes used the same event generator as the KTeV 
Kl — > e + e~77 measurement which included 0{a 2 ) 
radiative corrections in the calculation of the decay 
width [9(. Without radiative corrections the measured 
branching ratio would have been shifted by 1.7%. 

To implement 0(a 2 ) corrections, a cutoff of m 77 
1 MeV/c 2 was introduced to distinguish Kl — > e + e~7 
from the radiative decay, Kl — > e + e _ 77. For events 
below the cutoff, the second photon was not generated. 
The value of the cutoff was chosen so that the energy 
of the softer photon in Kl — > e + e _ 77 events near the 
cutoff would be too low to be detected. The detector 
acceptance for Kl — > e + e~7(7) events was found to be 
3.4%. For Kl — > tt !? ^ events, the acceptance was 
0.26%. 

For the purposes of the branching ratio measurement, 
the decay Kl — * e + e~7 can be defined as either all 
Kl — * e + e~7 and Kl — > e + e~77 events (the inclu- 
sive definition), or by introducing a cutoff to distinguish 
Kl — > e + e~7 from Kl — > e + e~77 (the exclusive defini- 
tion). We report our results using both of these defini- 
tions. For the exclusive result, Kl — > e + e~7 is defined as 
all events where the softer photon has an energy less than 
5 MeV in the kaon rest frame [E%£ < 5 MeV). This def- 
inition is chosen to coincide with previous measurements 
of the Kl — » e + e~77 branching ratio Q. 

The same definitional ambiguity also exists for the 
normalization mode. Since previous measurements of 
7T° — > e + e~7 have been inclusive [1 01 ] . using the inclusive 
definition as the normalization mode allows one to ex- 
tract more easily the absolute value of the Kl — * e + e~j 
branching ratio. 

Since the normalization mode had four more photons 



Branching Ratio Uncertainties 


% of BR 


A Can* 


A OlDIP 


Statistical Uncertainty 


0.33% 


0.030 


0.038 


BR(K L -» 7T u 7r u 7r^) Uncertainty 


2.83% 


N/A 


N/A 


Absolute 7 Inefficiency 


0.43% 


N/A 


N/A 


Drift Chamber Inefficiency 


0.37% 


0.009 


0.011 


Cut Variations 


0.33% 


0.013 


0.016 


Kaon Energy Spectrum 


0.23% 


0.011 


0.014 


Trigger Inefficiency 


0.21% 


N/A 


N/A 


Calorimeter Energy Resolution 


0.14% 


0.001 


0.001 


Background Level 


0.08% 


0.000 


0.000 


Detector Material 


0.07% 


0.008 


0.009 


Drift Chamber Hit Resolution 


0.04% 


0.002 


0.003 


0(a 3 ) Radiative Corrections 


0.03% 


0.008 


0.009 


Form Factor Dependence 


0.03% 


N/A 


N/A 


Total Systematic Uncertainty 


0.77% 


0.022 


0.028 



TABLE I: Uncertainties for the Kl —> e + e 7 branching ratio 
and form factor parameter measurements. 

than the signal mode, the largest systematic uncertainty 
was the absolute photon inefficiency in the calorimeter. 
The three sources of photon detection bias considered 
were the calorimeter geometry simulation, the simulated 
photon energy spectrum, and the electromagnetic shower 
containment. To measure the effect of the simulated 
calorimeter geometry, the outer edge of the calorime- 
ter was moved by 0.5 mm in the Kl — > 7r 7r 7r£, Monte 
Carlo events only. This resulted in a 0.226% variation 
in the branching ratio. The photon energy simulation 
was tested by shifting the Monte Carlo photon energies 
by 10 MeV, which caused the branching ratio to vary by 
0.219%. These variations were chosen based on detector 
survey data and information from other decay modes. 
Finally, an upper bound on the effect of imperfect elec- 
tromagnetic shower containment in the calorimeter was 
found by studying the low-end tail of the E/p distri- 
bution for electrons. Both Kl — > n e T v e events and 
Kl — > 7T 7r 7r° events were used to find best-fit shapes 
for the low-end E/p tail. Switching between these two 
shapes in the Monte Carlo caused a variation of 0.288% in 
the branching ratio. These three effects, added in quadra- 
ture, resulted in a 0.43% systematic uncertainty. 

The inefficiency of the drift chambers was also a source 
of systematic uncertainty. Two-dimensional inefficiency 
maps of each drift chamber were measured using Kl — > 
7r ± e =F ^ e decays. These inefficiency maps were then ad- 
justed by a constant factor to give the best agreement 
between data and Monte Carlo. One sigma variations 
about the best-fit value for this factor yielded a 0.37% 
systematic error on the branching ratio. 

An overall systematic uncertainty due to disagree- 
ments between the data and Monte Carlo was measured 
by studying the effect of varying the analysis cuts. The 
cut values on many important quantites were studied, 
such as the vertex position, E/p, energy in the photon ve- 
tos, reconstructed mass, pion probability, minimum pho- 
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ton energy, and total energy. These variations resulted 
in a 0.33% variation in the branching ratio. 

There was also a systematic uncertainty due to the 
simulation of the kaon energy spectrum. The ratio of the 
total kaon energy distributions between data and Monte 
Carlo exhibited a slope in both the signal and normal- 
ization modes. The slope was corrected by reweighting 
the Monte Carlo events based on total kaon energy. The 
reweighted Monte Carlo events were used for the central 
value measurement, and the variation in the branching 
ratio due to reweighting, 0.23%, was treated as a system- 
atic uncertainty. 

Several other smaller systematic effects have been eval- 
uated as well. The full list of systematic uncertainties is 
given in Table HI 

We have observed 83,262 Kl — > e + e^'j(j) events over 
a background of 79 events. Using 4,924,801 K L -> 
7r 7r°7r2) events to normalize the Kl — » e + e~7(7) rate, 
the inclusive {Kl —> e + e~7 + Kl — » e + e™77) and ex- 
clusive (E^2< 5 MeV) ratios of the two Kl decay widths 
have been measured to be 



T(e+e j)j 
r(7r 7r 7r0,) 
r(e + e~7) 



exc 



r(7r°7r°7r° ) 
r(e + e~7) 



trie 

exc 



T(Tr°Tr°TT Q D ) 



(1.3302 ± 0.0046 ± 0.0102) x 10~ 3 
(1.2521 ± 0.0044 ± 0.0097) x 10~ 3 
(1.2798 ± 0.0045 ± 0.0099) x 10~ 3 



exc 



Using the known values of BR(K L -> vr 7r 7r°), BR(tt -» 
77) and BR(7T° — * e + e~~f) inc [llj, the absolute Kl — * 
e + e _ 7 branching ratio is 

BR inc = (9.128 ± 0.032 ± 0.070 ± 0.252) x 10~ 6 
BR exc = (8.591 ± 0.030 ± 0.066 ± 0.238) x 10~ 6 

The third error listed is an external systematic error due 
to the uncertainty in the Kl — > 7r 7r 7r^, branching ratio. 
These results are in good agreement with all previous 
measurements with a factor of 3 improvement in the un- 
certainty on the ratio 



UaliaLij. 



Since both form factor models being considered are 
functions of the Dalitz variable x, the shape of the e + e~ 
mass distribution is very sensitive to the form factor 
parameters. Thus, each form factor parameter was ex- 
tracted by comparing the shape of the e + e~ mass spec- 
trum in data to several Monte Carlo samples with dif- 
fering form factor parameter values. A few such com- 
parisons are shown in Fig. [2] For each comparison, a 
bin-by-bin shape-x 2 value was calculated. A quadratic 
fit to these \ 2 values was used to determine each form 
factor parameter. This fitting procedure was performed 
separately for the two form factor parameters since the 
two models depend on x differently. 

The uncertainties associated with the form factor mea- 
surement are shown in Table |TJ Unlike the branching 
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FIG. 2: The variation in the shape of the e + e~ mass distribu- 
tion was used to determine the values of the form factor pa- 
rameters. Three comparisons between data (dots) and Monte 
Carlo samples (histograms) with different values of CaK* are 
shown as well as the final fit to the shape x 2 values 



s 
u 





} 










5 













E845 NA31 



NA48 KTeV KTeV 



KTeV 



1990 1990 1999 2001 



2002 2007 



FIG. 3: The graph shows a comparison of our measurement 
of Cax' with previous measurements. The outer error bars 
represent the total uncertainty and the inner error bars are 
the statistical uncertainty only (where applicable) [l2l , [l3l . 
15, 16] 



ratio uncertainties, the form factor uncertainties are re- 
ported as a variation in the measurement, not as a per- 
centage of the final result, since a measurement of zero 
has no special significance. The method for determin- 
ing these uncertainties was the same as for the branching 
ratio. 

The final fit for the form factor parameters Coik* and 
ctDiP yield 

Ca K * = -0.517 ±0.030 ±0.022 
aoip = -1.729 ±0.043 ±0.028 

The first uncertainty is from the % 2 minimization fit and 
the second is the total systematic uncertainty. 

Previous measurements of the BMS form factor have 
been reported in terms of ajp , where the constant pa- 
rameter C has been divided out. Since the expression for 
C involves several experimental quantities that have not 
been treated consistently in the past, we report CaK* 
to avoid confusion. Using previous measurements of ax* 
and their corresponding values of C, a comparison with 
our measured value of CaK* is shown in Fig. [3) In re- 
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cent years, the proper value for Colk- has been unclear 
due to a 3.1a difference between the two previous best 
measurements. Our measurement lies between these val- 
ues: 2.0c below the KTeV Kl — > ^i + /i~7 result and 2.4ct 
above the NA48 Kl — > e+e~7 result. 
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